Breakup of a nucleus with two weakly bound neutrons on a proton target: 

Single-scattering approximation 
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Three- and four-cluster breakup reactions in the ^^C- proton scattering are studied using three- 
body core-neutron-neutron model for ^®C. Single-scattering approximation (SSA) of four-particle 
equations for transition operators is used to calculate three- and four-cluster breakup amplitudes at 
200 and 300 MeV/nucleon energy near proton-neutron (pn) quasi free scattering (QFS) conditions. 
The differential cross section is sharply peaked at pn QFS point and decreases rapidly whenever 
kinematical conditions deviate from pn QFS. The accuracy of the SSA for the three-cluster breakup 
is estimated from a three-body model and is found to be as good as 6% at higher reaction energies 
in suitable angular regimes. Furthermore, under an additional approximation the three-cluster 
breakup amplitude factorizes into the pn transition operator and the overlap integral of two- and 
three-particle bound states. That approximation usually reduces the cross section, in some cases 
even up to 10%. 

PACS numbers: 24.10.-1,21. 45.-V, 25.10.-l-s, 25.60.Gc 



I. INTRODUCTION 

The structure and properties of exotic short-living 
nuclei are usually studied through their reactions with 
stable targets. On the theoretical side, the reliabil- 
ity of the extracted information depends on the valid- 
ity of the assumed interaction model and on the accu- 
racy of the theoretical methods used to solve the respec- 
tive few- or many-body problem. Among the standard 
tools for the analysis of breakup reactions are the dis- 
torted wave Born approximation, the eikonal approxima- 
tion [l| , time-dependent models Q , and the more sophis- 
ticated continuum-discretized coupled-channels (CDCC) 
method @. For reactions involving effective two-cluster 
systems like deuteron or one-neutron halo nuclei that are 
treated as three-body problems (two-cluster nucleus plus 
a target) also the exact Faddeev-type formalism 0, Q 
in the momentum-space framework has successfully been 
applied. Within this approach, once numerically well- 
converged results are obtained, all discrepancies with the 
experimental data can be attributed to the shortcomings 
of the used three-body model. However, practical appli- 
cations of the Faddeev-type formalism are limited so far 
to Hght projectiles and targets. 

On the other hand, Faddeev-type calculations can be 
used to check the accuracy of the traditional approximate 
nuclear reaction methods. Indeed, the comparison [f^ of 
Faddeev-type and CDCC results revealed that CDCC is 
a reliable method for deuteron-nucleus elastic scattering 
and breakup but fails for the breakup of a one-neutron 
halo nucleus, assumed to be a bound state of an inert 
nuclear core {A) and a neutron (n) , on a proton (p) tar- 
get. More specifically, the peaks of the differential cross 
section in energy and angular distributions of the core 
A are underpredicted by CDCC. Those peaks are due 
to pn quasi free scattering (QFS), i.e., only the neutron 
interacts with the proton and is knocked out from the 
{An) nucleus whereas the momentum and energy trans- 



fers to the core A vanish and, as a consequence, the core 
keeps the velocity of the incoming beam. Thus, near the 
pn QFS the reaction dynamics is dominated by the pn 
interaction which is very difficult to describe well when 
the three-particle CDCC wave function is expanded us- 
ing the nA bound and selected continuum states as done 
in standard CDCC (the transfer-to-continuum CDCC 
is more successful, but so far it is only available for three- 
body systems). 

For the same reason one may expect the failure of 
CDCC for the breakup of a core and two-neutron nu- 
cleus [Ann) on a proton target near pn QFS kinemat- 
ics, especially at higher energies relevant for the anal- 
ysis of present day and future experiments. Although 
four-body CDCC 0, Q and eikonal approximation [l[ 
calculations for breakup of two-neutron halo nuclei ex- 
ist, to the best of our knowledge, they were not applied 
in the pn QFS regime. Furthermore, additional techni- 
cal difficulties arise in CDCC if two-body bound states 
{An) exist On the other hand, the rigorous treat- 

ment of the four-particle scattering problem within the 
Faddeev-type framework was successfull y p erformed only 
for the four-nucleon system so far [Til. |T3|. The exten- 
sion of the Faddeev-type method to distinguishable par- 
ticles and high energies is very difficult technical chal- 
lenge, although in principle it may be feasible. On the 
other hand, such a huge effort may be unnecessary for 
the description of particular breakup reactions where the 
CDCC is expected to fail, i.e., for p -t- {Ann) breakup 
near pn QFS kinematics. The studies of ^^Be breakup 
on a proton target at 200 MeV/nucleon reaction energy 
around the pn QFS kinematics revealed that the 
single scattering approximation (SSA) reproduces quite 
accurately the results of the full Faddeev-type calcula- 
tions for not too small neutron scattering angles. Thus, 
one may expect the SSA to be reasonable also in the 
four-particle system at similar kinematical conditions. 
We therefore aim at developing the four-particle SSA to 
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study the breakup reactions p + {Ann) — )■ p + n + {An) 
and p + {Ann) -^p-\-n + n + A at higher energies. The 
numerical results use p + reactions as example. Un- 
fortunately, no experimental data is available. In the 
present work we will concentrate on theoretical aspects of 
the three-cluster breakup where one may create effective 
three-body model p + {Bn) —i'p + n + B with B = {An). 
Physically such a model is inadequate for weakly bound 
{An), but it allows to perform full three-body Faddeev- 
type calculations and, by comparing with the respective 
SSA, identify the kinematical regimes where the SSA is 
reasonable, and estimate its accuracy. 

In Sec. In] we derive the SSA for the breakup ampli- 
tudes. In Sec. [ni] we describe the employed interaction 
model. We study the three-cluster breakup in Sec. IIVI 
and four-cluster breakup in Sec. |Vl Summary is given in 
Sec. ED 

II. BREAKUP AMPLITUDES 

We consider the four-particle system interacting via 
short-range pairwise potentials vj where j labels the re- 
spective pair. We do not calculate the four-particle wave 
function explicitly but work in the momentum-space 
and use the integral form of the scattering equations 
as proposed by Alt, Grassberger and Sandhas (AGS) 
p^ . The AGS equations are equivalent to the Faddeev- 
Yakubovsky equations psj but are formulated for the 
four-particle transition operators Z^^p, i.e., 

U^p = {Go U Go)-^ ~5.p S,^+Y1 Go W;;, (1) 

The components of the operators are distinguished by the 
two-cluster partition and by the three-cluster partition; 
18 different combinations are possible. The two-cluster 
partitions, denoted by Greek letters, are either of 1-1-3 or 
2-1-2 type. The three-cluster partitions, denoted by Latin 
letters, are of 2-1-1-1-1 type and therefore are fully charac- 
terized by the pair of particles in the composite cluster. 
Obviously, the pairs i, j and k must be internal to the 
respective two-cluster partitions, i.e., i C p, j, A; C 7 and 
j C cr, whereas Sap = 1 — Sap- The components Z^^p with 
all allowed j and i describe the transition from the initial 
two-cluster partition p to the final two-cluster partition a 
. In our example of the p+ {Ann) scattering the initial 
two-cluster partition p is p + {Ann) with three internal 
pairs, i.e., {nn) and twice {An). Other two-cluster parti- 
tions of the 1-1-3 type are A + {nnp) and twice n + {npA) 
while those of the 2-|-2 type are {pA) + {nn) and twice 
{pn) + {An) . Note that the components of the transition 
operators exist even if the particles in the corresponding 
cluster do not bind as in the case of {pA) + {nn); those 
components contribute to breakup reactions. 

The free resolvent at the available system energy E is 
given by 

Go = {E + lO - Ho)-^ (2) 



where Hq is the kinetic energy operator. The two-particle 
transition operator 

tj — Vj + VjGotj (3) 

sums up the interactions for pair j up to all orders. Fur- 
thermore, all interactions within each two-cluster subsys- 
tem 7 lead to the respective subsystem transition opera- 
tors 

[/f - 1 5jk + U Go Ui^ . (4) 

i 

As derived in Ref. [13], the amplitude for the three- 
cluster breakup of the initial two-cluster state is 

($,|T,p|<fp) = ^(<i>,|C/f Gotfc GoU^^lW;). (5) 

The energy parameter in the operators of Eq. ^ is 
E = ep +Pp/2/ip with Cp being the energy of the bound 
state in the partition p and fip being the respective two- 
cluster reduced mass. The initial asymptotic channel 
state |<I>p) is a product of the bound state wave function 
in the partition p and the plane wave with momentum Pp 
between the two clusters. |$p) ~ J^i \4>p) is decomposed 
into its Faddeev components satisfying 

|0P = Go^%i,|0 (6) 

j 

and normalized such that ($o.|$p) = Sa-p 5{pa- — pp). The 
asymptotic three-cluster state \^j) is an eigenstate of the 
channel Hamiltonian Ho + Vj with the eigenvalue E. It 
is given by the bound state wave function for the pair 
j times two plane waves corresponding to the relative 
motion of three free clusters. 

The amplitude for the four-cluster breakup is taken 
over from Ref. [31, i.e., 

(<i>o|Top|<i>p) = Y.{^o\tj Go C/f Go h GoU^;\q^l). (7) 

jjki 

The four-cluster channel state |$o) is an eigenstate of Ho 
with eigenvalue E; it is a product of three plane waves 
(each is normalized to the Dirac i5-function) correspond- 
ing to the relative motion of four free particles. 

In the SSA, i.e., keeping only the terms of the first 
order in two-particle transition operators ([3]), the three- 
and four-cluster breakup amplitudes become 

(<i>,|7;^«^|<i>p)^^4p(<i>,|t.|<i>p), (8a) 

k 

($o|ToS;^|'i>p) = J2~^kp{<i>o\tk\<i>p). (8b) 

Here Skp is if /c C p and 1 otherwise. Thus, in Eqs. (|8]) 
the summation runs over all pairs that are external to 
the initial state partition p. Note that Eq. ([5a|) is the 
amplitude for the direct three-cluster breakup and not 
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for rearrangement breakup, i.e., the final bound pair j is 
internal to the initial partition p. 

In the following we consider the initial partition p to be 
of the 1+3 type, i.e., 1(234) with particle 1 as spectator 
and the bound state of particles (234). In such case the 
amplitudes (|8|) have three contributions J2k ^kptk — ti2 + 
^13+^14 corresponding to the interactions of the spectator 
particle 1 with each particle in the cluster (234) but no 
interactions within the cluster. 

We start with the four-cluster breakup amplitude (jSbp . 
Let 771(1 be the mass of the particle a and its final-state 
momentum with a — 1,2,3,4. Furthermore, the initial 
momentum of particle 1 and of cluster (234) we denote 
by and k^, respectively. Obviously, k^, k^ and ka 
are related by momentum and energy conservation, i.e.. 



E ka = kl + kl = K and E kl/2ma = 



TO2 + + m4. 



k{ /2mi + 
We do not 



hp /2mp — E with 
assume a particular frame, thus, the results are valid with 
the energy E and momenta given in any frame. The 
explicit momentum-dependence of the ti2 term is 



($0|il2|*p) - (pi2Kl2(ei2+zO)|p'i2)(q2P34|*p) (9) 



with 



Pab 



Pl6 

kl 



^2 = 



mfaka ~ mgkb 

771a + mb 
rubkl - TOikfc 

rrii + nib 
ki +k6-ki, 

— [{ms + 7714 )k2 - "12 (ks + k4)] 

771„ 



(10a) 

(10b) 
(10c) 
(lOd) 



Alternatively, another set of Jacobi momenta for the 
(234) subsystem could be chosen to represent |$p). Here 
and in the following the wave functions with notation- 
ally indicated momentum dependence refer to the inter- 
nal motion of the respective bound cluster; the part cor- 
responding to the free motion between clusters is taken 
out. The pair transition operator ti2 depends on the en- 
ergy available for the relative motion of particles 1 and 
2, i.e., 



ei2 = E - 



(ki+k2)^ 

2m3 2m4 2(7711+7712) 



^4 



(11) 



Due to momentum and energy conservation ei2 — 
^^2/2^12 with the reduced mass Pab = Wa"ifc/('7^a + "*■&); 
thus, the pair transition operator ti2 has to be evalu- 
ated half-shell. Furthermore, although for brevity not 
explicitly indicated in our notation, the breakup ampli- 
tudes ^ and pair transition matrices (jSj are operators 
in the spin-space implying summations over all interme- 
diate spin states. In Eq. Q this summation involves only 
the initial spin projection of the particle 2. 

The momentum-dependence of the ti3 and ti4 terms 
has the same structure and can be easily obtained from 
Eq. ^ by the respective permutation of particle labels. 



Considering the three-cluster breakup we assume that 
particles 3 and 4 build the final bound pair with mass 
771b = 7113 + "^4, total momentum ks, internal wave 
function (p34|<I's) and energy e^. The kinematics of the 
QFS condition for the pair of particles 1 and 2 reads 
ks « k§^^ = {mB /iTipj'kp. The most important con- 
tribution of the breakup amplitude is the <i2 term rep- 
resented diagrammatically in Fig. [T] (a). It includes the 
interaction between the spectator and the struck particle. 



(*B|il2|*p)= J rf^P34[(Pl2|il2(ei2 +«0)|p'i2) 
X (*B|P34)(q2P34|*p)]- 



(12) 



The definitions of the involved momenta as given in 
Eqs. (fTOl) are valid also here but k3 = 7n3kB/(7n3 + 7)i4) + 
P34 and k4 = r7i4kB/(77i3-|-77i4) — P34 are not independent 
variables anymore. The momentum and energy conser- 
vation changes to ki + k2 + k^ = k^^ + 

kj/2mi + fc|/2ni2 -I- kg/2mB = fp 
p = E. The energy ei2 as defined in Eq. pT|) is 



k\'/2m^ 



k;, = K and 

+ fcf /27711 -I- 



valid as well but it is more appropriate to express it via 



Pah, I.e., 



, _ P?2 , ^ P34 
ei2 - 7; h es - T, • 

Api2 Z^34 



(13) 



Thus, this time the two-particle transition operator has 
to be evaluated fully off-shell. More importantly, Eq. 
demonstrates that ^12(612 -I- iO) depends on the integra- 
tion variable P34 and therefore cannot be taken out of the 
integral in Eq. ((T2|) . An additional approximation in the 
energy-dependence of ti2 is needed to factorize the three- 
cluster breakup amplitude ([T2]) into ti2 and the overlap 
integral (01) 



= / d^P34{^B\p34){ci2Pu\^p 



(14) 



We introduce the overlap integral approximation (OIA) 
of the three-cluster breakup amplitude as 



OIA 
Bp 



(Pl2|tl 



P?2 



2p 



12 



*0)|p'i2>X^(q2) 



(15) 

Under this particular approximation the two-particle 
transition operator needs to be evaluated half-shell only 
and the four-particle SSA becomes formally a three- 
particle SSA since the amplitude (fT5|) has exactly the 
structure of the SSA breakup amplitude for the three- 
particle system (l-l-2-l-B) where the (2B) bound state 
wave function is replaced by the 01 Xp (q2)- 

The ^13 and tu terms of the three-cluster breakup am- 
plitude (|8a|) have the structure diagrammatically repre- 
sented in Fig. [T] (b). They include the interactions be- 
tween the spectator and the particles that remain bound 
in the pair B, i.e.. 



i^Blhsl^p) = J d3p34[(pi3|il3(ei3 +«0)|p'i 
X (*i3|P34)(q2P34|«'p)] 



(16) 
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(a) (b) 

FIG. 1. (Color online) Two types of contributions to the 
three-cluster breakup amplitude in the SSA. The two-particle 
transition operator is represented by a box while two- and 
three-particle bound states are represented by filled arcs. The 
diagrams (a) and (b) correspond to Eqs. (|12|l and (I16|) . respec- 
tively. 

with q2 = [(ms + m4)k2 - m2(k3 + ki)]/mp, P34 = 
(m4k3 — m3k4)/(m3 + 7124), and 643 defined according 
to Eq. pT|) with the permutation (2 o 3). Thus, all 
relative momenta in Eq. ()16p depend on the integration 
variable P34 and therefore no simplifications are possible. 
However, our numerical calculations revealed that in the 
regime relevant for SSA, i.e., near pn QFS, the contri- 
butions of the type are very small as compared to 
and can be safely neglected. There are good physics 
reasons for this since, as mentioned, tis and ti4 terms 
describe the spectator interaction with the particles that 
form the bound pair B but not with the struck parti- 
cle. The suppression of contribution (I16p comes from 
high momentum components of the involved wave func- 
tions while in Eq. ()12[) low-momentum components are 
decisive. 

For p + (Ann) breakup reactions considered in this 
work two additional remarks regarding the symmetriza- 
tion and the long-range Coulomb interaction are needed. 

First, the channel states in Eqs. ([8|) have to be anti- 
symmetric under exchange of the two neutrons. For the 
initial state |$p) this is achieved by including only an- 
tisymmetric two-neutron partial waves and keeping only 
two independent Faddeev components j^p when solving 
the three-particle bound state problem. The antisym- 
metrized final asymptotic three- and four-cluster states 
to be used in Eqs. ^ are 

\K) = ^il~Pnn)\'f.) (17) 

with Pnn being the two-neutron permutation operator. 
Since |<i>p) is already antisymmetrized and ^12+^13+^14 is 
symmetric under exchange of the neutrons, the antisym- 
metrization (jl7p simply yields an additional \/2 factor 
for the amplitudes ([5]). 

Second, the AGS equations for the transition opera- 
tors are defined only with short-range potentials Vj . The 
Coulomb interaction (with no more than two charged 



clusters) can be included using the method of screening 
and renormalization [T9I-I2H but only in the full form of 
AGS equations where the unscreened limit for the renor- 
malized amplitudes exists. This is not the case in the 
SSA: tpn being Coulomb-free needs no renormalization 
while tpA includes Coulomb and needs renormalization 
half-shell To']. Thus, strictly speaking, their sum has no 
unscreened limit neither with nor without renormaliza- 
tion. Practically, this is not a problem since, according 
to our calculations, near the QFS conditions where the 
SSA is expected to be reasonable, the contribution of tpA 
is very small and can safely be neglected. 



III. INTERACTIONS 

We takep+^'^C^ p + n + ^^G andp+^^C^ p + n + n + 
^^C reactions in inverse kinematics as a working example 
in the numerical calculations of this paper. The ground 
state of ""^^C, the core A, is well separated from the excited 
states (6.093 MeV) while ^^C is weakly bound (es = 
— 1.218 MeV) one-neutron halo nucleus for which a simple 
two-body model of core A and neutron n is assumed to 
be adequate. The neutron separation energy of ^^C is 
4.250 MeV such that core and two-neutron model of ^^C 
ground state 0+ with = —5.468 MeV appears to be 
quite reasonable. 

As the pn and nn interactions we take the charge- 
dependent Bonn (CD Bonn) potential [l^. The nA 
interaction is taken over from Ref. [2^; it has central 
and spin-orbit parts of Woods-Saxon shape adjusted to 
^^C bound states 2si/2 and 1^5/2 and to ^^C neutron 
separation energy in lpi/2. The deeply-bound Pauli 
forbidden states lsi/27 lPi/2; a-nd IP3/2 are projected 
out as described in Ref. [23i]. The optical potential as 
parametrized by Koning & Delaroche ^24] plus screened 
Coulomb is used for pA interaction; one could probably 
find a better parametrization but, as already mentioned, 
this is irrelevant since tpA yields negligible contribution 
to SSA near pn QFS. 

The above nn and nA potentials alone do not yield 
the desired value €p = —5.468 MeV for the ^^C ground 
state. A simple way to remedy this shortcoming is to add 
a three-body force (3BF) acting in the nnA subsystem 
only. In this way it does not affect the functional form of 
the SSA breakup amplitudes ([5]). The coordinate-space 
calculations usually take a central 3BF depending on hy- 
perradius. Our calculations are in momentum space so 
we take a central 3BF depending on hypermomentum /C 
that in the three-particle (234) subsystem has the form 

(k2k3k4|iymk^k;) = (4^)2u,3.g(/C2).g(/C'') (18) 

with /C^ — mAr[A:|/m2 + k^/m^ + k\/m4 — (k2 + k3 + 
k4)^/mp] and JC' defined analogously where mjv is the 
average nucleon mass. Note that JC^ /2mM is the internal 
motion kinetic energy of the (234) subsystem. We chose 
the form factor gilC"^) = exp (— /C^/2A^) as Gaussian. 
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The form ([T8)) of the 3BF is inspired by the effective field 
theory . Once for each value of the cutoff parameter 
A the strength W3 is adjusted to reproduce Cp — —5.468 
MeV, the predictions become practically independent of 
A. For example, changing A by a factor of 2, from 2 fm^^ 
to 4fm^^, yields less than 3% changes in the breakup 
cross sections. Our standard choice is A = 3fm~^. 

With 3BF included the Faddeev components of the 
three-particle bound state obey the equation 

Wp)^GoJ2s,,t,\<lyl)+Goil + kGoHwY,\<l>l) (19) 
j 3 

with 'Y^rji = 1. Different rji choices correspond to differ- 
ent splitting of the 3BF contributions among the Faddeev 
components but yield identical ep and |$p). 

We calculate the two-particle transition operators tj 
and the bound state wave functions |$_b) and \^p) in the 
momentum-space partial-wave basis but then transform 
them into the plane-wave representation as needed in 
Eqs. ([MTBl). To obtain converged results the two-particle 
interactions Vj have to be included up to the two-particle 
relative orbital angular momentum i^^ax- We find that 
^max = 11 for pn, 3 for nn, and 3 for nA is sufficient. 
The test calculations proving the negligible contribution 
of tpA used i^iax = 20. 

We note that at ei2 = ed with Cd = -2.223 MeV be- 
ing the CD Bonn prediction for the deuteron bound state 
energy, the tpn transition operator in the ^5'i — par- 
tial wave exhibits the deuteron bound state pole. In the 
integrals it is treated by the subtraction technique. 

IV. THREE-CLUSTER BREAKUP 

We consider the three-cluster breakup where two clus- 
ters (a and b) are detected with momenta and k;, (all 
single-cluster momenta in this section refer to the lab 
frame). The momentum of the undetected cluster c is 
fully determined by the momentum conservation. The 
energy conservation renders ka and kb not independent; 
for a fixed ka there may be up to two solutions for k^, and 
vice versa. The five independent kinematic variables for 
the fully exclusive fivefold differential cross section are of- 
ten chosen as the polar and azimuthal scattering angles 
= (0a,'/'a) ^nd = (©6, "ySb) of the two detected 
particles and one energy i?a, i.e., 

d^cr {2Tr)^ma'm},mckak'^ 
diladflbdEa V\{mb + mc)kl - mf,(K - ka) • kb| 

xl5]|($s|rsp|$p)p. 



Here V = jk^^/mi — k^/mpj is the incoming flux. The 
sum runs over all initial and final spin states, while gi — 
(2si-|-l)(2sp-|-l) takes care of the spin averaging in initial 
state, si (sp) being the spin of the particle 1 (cluster p). 



In our example of p -I- ^^C— > p + n-\- C reaction we as- 
sume that beam of ^^C is impinging on target p and the 
detected particles are ^^C in its ground state {B) and p. 
We compare the differential cross section (pUj) calculated 
in four different ways, depending on the scattering am- 
plitude ($_b|Tbp|$p). The four-particle SSA as given by 
Eqs. (|12ill3p is labeled SSA-4b. Its approximation (US]) is 
labeled 0IA-4b. The importance of higher order interac- 
tions between the three clusters can be estimated by cre- 
ating an effective three-body model where the composite 
cluster B is treated as a single inert particle. The nB po- 
tential is real in 0"*" wave and supports bound state with 
energy Cp — cb, while t he p B interaction includes opti- 
cal Koning & Delaroche [2J] and Coulomb potentials; the 
latter is treated using the method of screening and renor- 
malization [23| . The results obtained by solving full 
three-body Faddeev- type equations, i.e., formally sum- 
ming up multiple scattering (MS) series up to infinite 
order, are labeled MS-3b. The SSA of this model in- 
cluding only tpn term is labeled SSA-3b. Although such 
model makes physically little sense owing to halo nature 
ofB, but the ratio [(MS-3b)-(SSA-3b)]/(MS-3b) should 
be a reasonable accuracy estimate of SSA-4b. Note that 
for comparison the results of the three-body model are 
multiplied by 2 to account for the two neutrons in the 
original four-particle model. 

We concentrate on the kinematic regime near pn QFS, 

i.e., \s.B ~ — {vtlb /'mp)\i^p. In terms of angles 

ops 

and energy this means Qb ~ and Eb ~ E]^ — 
{k'^^Y /2mB = {mB/mp)Ep where Ep is the beam en- 
ergy of ^^C. Note that at these conditions also q2 vanishes 
in Eq. 

In Fig.[2]we show the results at £^p/16 = 300 MeV. We 
fix Qb = 0°, 'fB = 0°, ipp = 180°, and vary Op. In all 
used approaches the differential cross section peaks quite 
sharply at (or very near to) Eb = £'2^^. This is due 
to sharp increase of the s-wave components of the bound 
state wave functions for vanishing relative momenta. At 
Op = 15° and 30° there is significant difference between 
MS-3b and SSA-3b indicating that the SSA is not reli- 
able in this regime. On the contrary, at Op = 45° and 
60° the agreement between MS-3b and SSA-3b gets bet- 
ter, of the order of 6% at the peak. Thus, for these 
kinematic configurations SSA-4b should be of a compa- 
rable accuracy. The SSA-4b results, taking into account 
the three-particle structure of ^^C, are considerably lower 
than the ones of SSA-3b where this aspect is neglected. 
The 0IA-4b, involving an additional approximation in 
the energy-dependence of underestimates the SSA- 
4b up to 9%, most sizably at Op = 45°. 

In Fig. [3] we fix Op = 45° where the SSA is expected 
to be at its best but vary Qb and include also results 
for -Bp/16 = 200 MeV. At this lower energy the dif- 
ference between MS-3b and SSA-3b is more significant, 
around 15% at the peak while at E^/IQ = 300 MeV it 
remains below 10%. Thus, these results confirm once 
more that SSA is more reliable at higher energies. On 
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FIG. 2. (Color online) Differential cross section for 

^^C(p,pn)^^C reaction in inverse kinematics as a function of 
the final ^^C energy Eb- The energy of the beam is 
300 MeV/nucleon. ^^C and proton are detected at angles 
(Ob = 0°,v5s = 0°) and (Qp = 15, 30, 45, 60°, (^p = 180°). 
Resuhs of SSA-3b (dotted), MS-3b (dotted-dashed), SSA-4b 
(solid) and 0IA-4b (dashed) calculations are compared. 
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the other hand, sensitivity to the treatment of ip„ energy- 
dependence is slightly more pronounced at lower energies: 
the difference between SSA-4b and 0IA-4b reaches 11% 
at E^'pjX'o = 200 MeV as compared to 9% at £;^/16 = 300 
MeV. A further message from Fig. [3] is that at higher en- 
ergy the differential cross section is more sharply peaked 
around the QFS point as it decreases faster with increas- 
ing Ob- 

In Fig. H we show the results at -B*/16 = 300 MeV 
where a neutron is detected instead of a proton. As in 
Fig.Hwefixes = 0°, = 0°, = 180°, and vary e„. 
In this case the best agreement between MS-3b and SSA- 
3b, about 6%, is observed at intermediate angles 6„ = 
30° and 45° while at e„ = 15° and 60° the difference 
gets above 20%. The effect of OIA is most sizable at 
e„ = 45° reaching 8%. 

As can be concluded from the above results, the agree- 
ment between MS-3b and SSA-3b, i.e., the reliability of 
the SSA, depends not only on the reaction energy but 
also on the kinematical configuration. This dependence 
can be understood by inspecting the nB and "pB relative 
energies in the final three-cluster breakup state, EnB and 
E^B ■ The values corresponding to QFS peaks in Fig. [2] 
are E^b = 258, 207, 137, and 66 MeV and E^b = 19, 
70, 141, and 211 MeV for Op = 15, 30, 45, and 60°, re- 
spectively. The values for Fig. 2] are obtained by simply 
interchanging E^b and E^b- Thus, the SSA appears to 
fail when either E^b or E^b is not large enough, with 
EpB being more decisive, possibly due to -pB Coulomb 
interaction. A detailed investigation of the MS-3b con- 



FIG. 3. (Color online) Same as Fig. Hbut for ^®C beam en- 
ergy of 200 (left) and 300 (right) MeV/nucleon and scattering 
angles (Gs = 0.0, 0.2, 0.4°, <^s = 0°) and (Sp = 45°, (^p = 
180°). 
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tributions reveals that, depending on EpB and EnB, its 
difference to SSA-3b is dominated by 2nd order terms 

tpsGotpn or tnBGotpn- 

We also note that at the QFS peaks the relative pn 
energy is around 150 MeV; thus, it is still in the regime 
where the underlying pn potential is well constrained by 
the experimental data. 

The OIA as given by Eq. (fT5|) results in higher values 
of the energy ei2 at which the pn transition operator tpn 
is evaluated as compared to the original SSA (IT^ . This 
might lead to smaller tp„ values in OIA since the matrix 
elements of ip„ in average decrease with increasing en- 
ergy much like the total pn cross section does. This may 
explain why the 0IA-4b in most cases underestimates the 
SSA-4b predictions. Furthermore, we conjecture that the 
use of the overlap integrals in many-body reaction models 
would have qualitatively similar effect, i.e., would reduce 
the predicted cross sections. 
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V. FOUR-CLUSTER BREAKUP 

A kinematically complete measurement of the four- 
cluster breakup requires the detection of three-clusters at 
least. As this is extremely difhcult, semi-inclusive observ- 
ables like momentum distributions are usually measured. 
Nevertheless, here we present results for the fully exclu- 
sive differential cross section that serves as a intermediate 
step for the calculation of semi-inclusive cross sections. 
With this goal in mind it is advantageous to choose the 
relative momenta as kinematical variables, i.e., 



nip + 2mn 



(21a) 
(21b) 



_ (rrip -I- 2TO„)kA - mA(k 

where the subscripts ni and 712 distinguish between the 
two neutrons and M = wla + nip + 2m„. The associ- 
ated relative energies are E,^, = fc^/2/^a;, Ey — ky/2fiy, 
and Ez = kl/2fiz with reduced masses = nin/2, 
pLy = 2ninnip / (nip + 2nin) , and = niA(nip + 2nin)/M. 
For example, kz and Ez^z/niA are the momentum and 
energy of the nuclear core A in the four-particle center- 
of mass (cm.) system, while k^^ and E^ are the two- 
neutron relative momentum and energy. Owing to mo- 
mentum and energy conservation there are eight inde- 
pendent kinematic variables; we choose them as the po- 
lar and azimuthal scattering angles ^Ij = {Qj,ipj) with 
J = x,y^z and two energies E^ and Ez. In this rep- 
resentation the spin-averaged eightfold differential cross 
section is 



(2^)' V- 



dilxdilydflzdErcdEz 



Vgi ^ 

X fJ'xtJ'yfJ'zkxkyk 



(22) 



FIG. 5. (Color online) Differential cross section for 

^^Cip,pnny'^C reaction at El/16 = 300 MeV as a func- 
tion of the relative nn energy Ex for selected values of Ez 
near pn QFS kinematics and angles 6^;, Oy, and Oz whereas 

ifx = tPy = Ifiz = 0°. 



where, as in Eq. (j20p . the sum runs over all initial and 
final spin states. The single-particle cross section for the 
core A in the cm. frame can simply be obtained as 



dVlAdEA 



niA 



d^a 



d^xd^lydEx 

fj,z dilxditydilzdExdEz 



The pn QFS implies k^ 



.QFS 



(23) 

(TO^/mp)k* and 



k„2 fa ^n2^ — (™n/™p)kp but also vanishing relative 
momentum Pn2A = {mA^7i2 ~ 'n^n'^A)/{niA + nin) and 
relative energy En2A (here for simplicity we assume that 
ni is knocked out, but in practical calculations the am- 
plitudes are antisymmetric with respect to the two neu- 
trons). Thus, due to low relative nA energy the SSA is 
expected to be less accurate than for the three-cluster 
breakup in suitable kinematics. Nevertheless, in Fig. [5] 
we show the results for the fully exclusive differential 
cross section ([22| of the p + ^^C— ?> p + n + n + "^^C reac- 
tion at E^p/IQ = 300 MeV. We fix the azimuthal angles 
ifx = fy — = 0° and vary the polar angles near 
a pronounced QFS peak at Qx — 35°, Qy — 90°, and 
Qz =0°. The cross section is shown as a function of the 
relative nn energy Ex for selected values of Ez around 
£;QFS ^ [mlmy [nznipM^)]El. The resuhs demon- 
strate that the differential cross section decreases rapidly 
whenever kinematical conditions deviate from QFS. 

Finally, we note that the reliability of the SSA pre- 
dictions for the semi-inclusive cross section (|23| may be 
even more limited. The integration in Eq. (|23p unavoid- 
ably includes regions of phase space with low relative pA 
energy where the higher order terms omitted in the SSA 



8 



may be significant. 



VI. SUMMARY 

Starting with full four-particle AGS equations for the 
transition operators we derived three- and four-cluster 
breakup amplitudes in the single-scattering approxima- 
tion. Numerical calculations were performed for ^^C 
breakup on a proton target with the core and two- 
neutron model for the ^^C nucleus. Breakup reactions 
at 200 and 300 MeV/nucleon energy near pn QFS con- 
ditions were studied. In the case of the three-cluster 
breakup an additional thrcc-body model, although being 
physically inadequate, allowed to estimate the accuracy 
of the SSA and thereby identify the kinematical regimes 
where the SSA is reliable. These regimes correspond to 
higher reaction energies and to higher relative energies 
between the composite cluster and any of the nucleons 
and are realized at proton or neutron scattering angles 
around 45°. There the accuracy of the SSA becomes as 
good as 6%. Furthermore, we have shown that an addi- 
tional approximation in the energy dependence of the pn 
transition operator is needed to factorize the SSA of the 
three-cluster breakup amplitude into tpn and the overlap 



integral of two- and three-particle bound states. This 
approximation usually reduces the cross section, in some 
cases even up to 10%. 

No SSA reliability test was possible for the four-cluster 
breakup but, given the conclusions drawn in the three- 
cluster case, it is expected to be less accurate. An exten- 
sion of the present calculations to include also the double- 
scattering terms, especially those between the core and 
neutron at low relative energy, would be highly desirable. 

The obtained three- and four-cluster breakup results 
demonstrate that the differential cross section is sharply 
peaked at pn QFS point and decreases rapidly whenever 
kinematical conditions deviate from pn QFS. 

The present numerical calculations are limited to p -|- 
^^C reactions. However, the formalism is applicable also 
to other nuclei like ^^Be, ^°C, and ^^O, and, in the case 
of the four-cluster breakup, also to ^He, ^^Li, and ^^C. 
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